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Abstract 

We study the effects of oblateness up to J4 of the primaries and power-law density profile 
(PDP) on the linear stability of libration location of an infinitesimal mass within the frame¬ 
work of restricted three body problem (R3BP), by using a more realistic model in which a 
disc with PDP is rotating around the common center of the system mass with perturbed 
mean motion. The existence and stability of triangular equilibrium points have been ex¬ 
plored. It has been shown that triangular equilibrium points are stable for 0 < /i < /Xc 
and unstable for /Xc < /U < 1/2, where Hc denotes the critical mass parameter. We find 
that, the oblateness up to J 2 of the primaries and the radiation reduces the stability range 
while the oblateness up to J4 of the primaries increases the size of stability both in the 
context where PDP is considered and ignored. The PDP has an effect of about ~ 0.01 
reduction on the application of /Xc to Earth-Moon and Jupiter-Moons systems. We find 
that the comprehensive effects of the perturbations have a stabilizing proclivity. However, 
the oblateness up to J 2 of the primaries and the radiation of the primaries have tendency 
for instability, while coefficients up to J4 of the primaries have stability predisposition. In 
the limiting case c = 0, and also by setting appropriate parameter(s) to zero, our results are 
in excellent agreement with the ones obtained previously. Libration points play a very im¬ 
portant role in space mission and as a consequence, our results have a practical application 
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in space dynamics and related areas. The model may be applied to study the navigation 
and stationkeeping operations of spacecraft (infinitesimal mass) around the Jupiter (more 
massive) -Callisto (less massive) system, where PDF accounts for the circumsolar ring of 
asteroidal dust, which has a cloud of dust permanently in its wake. 

Keywords: Restricted three-body problem; Libration points; Oblateness; Power-law density 
profile. 


1 Introduction 


Three-body problem is an important class of problem in classical and quantum mechanics which 
involves modeling the motion of three particles or massive bodies subject to their mutually 
perturbing gravitational attractions, e.g., the Sun-Earth-Moon system. For this system, the 
Sun’s mass is so dominant that it can be treated as a hxed object and the Earth-Moon system 
is treated as a two-body system from the point of view of a reference frame orbiting the Sun 
with that system. The simplest form of th e three-body problem is the restricted three-body 
problem (R3BP) fISingh and Taura I . l2014al) . It describes the motion of an inhnitesimal mass 


moving under the gravitational influence of two massive bodies known as the primaries which 
move in a circular orbits around their center of mass on the account of their mutual attractions. 
Moreover, the inhnitesimal mass is not inhuenced by the motion of the primaries. 

To tackle the problem of classical restricted three-body problem, Lagrange considered the 
behavior of the distances between the bodies without Ending a general solution. He obtained 
numerous equations from which he discovered two distinct classes of constant-pattern solutions. 
The first is referred to as the collinear (in which one of the distances is the sum of the other two) 
and the second is equiangular (in which the three distances are equal). The two distinct classes 
result into Lagrangian points Li, L2, .^3, L4 and L5, respectively which form the hve special points 
in rotating frame of reference where gravitational equilibrium can be maintained. 

The triangular libration points of the dynamical system denoted by L4 and L5 are sometimes 
called as the triangular Lagrange points or Trojan points. The name “Trojan points” comes from 
the Jupiter Trojans at the Sun-Jupiter L4 and L5 points, named after characters from Homer’s 
Iliad. The Sun’s gravity perturbs objects in the libration points so as not to be a perfect picture, 
but rather fairly close. The distances to mi and m2 are equal at L4^5 and thus, make the points 
balance. The resultant force acts through the barycenter of the system due to the gravitational 
forces from the two massive bodies which are in the same ratio as the masses of the primaries. 

In Figure [H we show the position of these Lagrangian points in a Sun-Earth system. These 
hve Lagrangian points have very important applications in astronautical. A great number of 
space missions have been completed while some operations are still in progress. The Solar and 
Heliospheric Observatory (SOHO) lunched in 1995 and Microwave Anisotropy Probe (MAP) 
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lunched in 2001 by NASA are currently in operation Sun-Earth Li and L 2 respectively. Solar 
TErrestrial RElations Observatory-Ahead (STEREO-A) made its closest pass to L 5 recently, 
on its orbit around the Sun. Asteroid 2010 S016, is currently proximal to L 5 but at a high 
inclination. In view of the importance of Lagrange points to the exploration and development 


of space, the dy namics an c 


influence in ref. f Trivailo 


stab ility of a satellite were studied under multiple Trojan asteroids 


200711 . 



Figure 1: (Color online) Position of Lagrangian points in a Sun-Earth system 


In studying Classical R3BP, it is required to put into consideration the oblateness of a body 
because, some bodies in our solar system for instance the earth and some in stellar system such 
as Altair are not spherical but oblate. In light of this, an enormous number of studies have been 
made in the recent years. The literature is huge, however, the following references give some 
older as well as recent studies. Letelier and Silva obtained a particular solution to the restricted 
three-body problem where the bodies ar e allowed to either lose or gain mass to or from a static 
atmosphere ( Letelier and da Silva I. 2011 1 . The motion of an inhnitesim al body in the generalize d 
restricted three-body problem have been presented by Singh and Taura (jSingh and Tanra 1.120131) . 
It was found that, in addition to the usual hve equilibrium points, there appear two new collinear 
points Lni, Ln 2 due to the potential from the belt, and in the presence of all these perturbations, 
the equilibrium points Li, L 3 come nearer to the primaries; while L 2 , L 4 , L 5 , Lni move towards 
the less massive primary and Ln 2 moves away from it. 

The dynamics of galactic systems with central binary black holes was studied by modify¬ 


ing the re stricted three body pr oblem, in which a galactic po tential was added as an ex 


potential fl.Iiang and Yeh 1201411 . Narayan and Shrivastava flNaravan and Shrivastava 


iernal 


201311 
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discussed the effects of oblateness and the photo-gravitational of the primaries on the location 
and the stability of the triang ular equilibrium points in the ell iptical restricted three-body prob¬ 
lem. Abouelmagd and Sharaf f Abonelmagd and Sharaf . 20131) explored the existence of libration 
points and their linear stability when the more massive primary is radiating and the less massive 
primary is an oblate spheroid. The existence and linear stability of equilibrium points in the 
perturbed Robe’s circular restricted three-body problem under the assumption that the hydro- 
static equilibrium figure of the hrst primary is an oblate spheroid has been examined recently 


(ISingh and Cvril-Okeme 


,| 2 ^. 


Analytical study of the dynamics of the third body within the framework of R3PB hav¬ 
ing a variable mass which changes accordingly to Jeansa la w has bee n presented r ecent ly by 


Abouelmagd and Mostafa (lAbouelmagd and Mostafa 1 . 120151) . In ref. (iKholostova 


20151). the 


author investigated the stability of triangular libration points in a planar restricted elliptic three- 
body problem for two sets of parameters corresponding to the cases of double resonances. The 
existence of the libration points and their linear stability within the framework of the R3PB, 
considering the effects of the hrst two even zonal harmon i cs pa rameters with respect to both 


primaries have been studied recently (lAbouelmagd et ah l2015f) . Couple R3B P is one of the 


examples of restricted four-body problem (R4BP). Within this context, Falaye flFalave 120151) 
studied the motion of an inhnitesimal mass by assuming that the primaries of the system are 
radiating-oblate spheroids surrounded by a circular cluster of material pointS; _ 


In this research work, we intend to extend the work of Singh and Taura fISingh and Taura 


2014al) by considering a modihed restricted three body problem without the effect from circular 
cluster of material point but with a disk, which rotates around the cente r of mass of the systeni 


with perturbed mean motion and also, the work of Kushvah and Kishor flKushvah et ah 120121) 


by examine the effect of oblateness of the primaries up to J 4 of the zon al harmonics on the stabilit y 


of equilibrium points. Our results will also generalize the recent work (lAbouelmagd et ah 


20151 


2 Mathematical formulation of the model 


Let rrii {i = 1,2) denote the masses of the primaries (mi for the more massive primary and m 2 for 
the less massive primary) and let the mass of an inhnitesimal body moving in the plane of motion 
of the primaries be m. The positions of the primaries are dehned relative to a rotating coordinate 
frame oxyz whose x-axis coincides with the line joining mi and m 2 and whose origin coincides 
with the center of mass of mi and m 2 ; while y-axis is perpendicular to the x-axis, the 2 ;-axis 
is perpendicular to the orbital plane of the primaries. Let ri and r 2 be the distances between 
m and the primaries mi and m 2 respectively, and R be the distance between the primaries, m 
moves in orbital plane under their mutual gravitational helds. The sum of m i and m^ is unity and 


h = ’^ 2 / (mi -|-m 2 ) is the mass ratio with the distance between taken as one (lAbouelmagd et al 
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2ni4J) . Furthermore, the unit of time is chosen to make both the constant of gravitation and 
the unperturbed mean motion equal to 1. To study the position of the inhnitesimal mass in the 
plane of motion of the primaries, we can either use the sidereal system of coordinates, or the 
synodical system of coordinates. In a synodical system of coordinates, we write the coordinates 
of mi, 1712 and m as (—/i, 0, 0), (1 — /i, 0, 0) and {x, y, z) respectively. Hence, the kinetic energy 
(K.E.) of the inhnitesimal mass in this barycentric coordinate system rotating abont z-axis with 
constant angular velocity n, is represented as 


K.E. = -m + y'^) + 2n{xy — yx) + (x^ + y‘^)\ , 


( 1 ) 


where the dot denotes differential with respect to time. Now, we take the potential energy of the 
inhnitesimal mass under radiating-oblate primaries as 

3H2\ fl Bi 3 B 2 ^ 


V = -Gm 


9i(i - f*) I 


1 Hi 

+ 1 — + 


2 rf 


Srf 


( 2 ) 


Furthermore, we consider a more realistic model in which a disk with PDF, is rotating around 


the common center of mass of the system. The PDF of the disc having a thickness h 10 ^ is 


p{r) = cr P f Kushvah et ah . 12012 1. where p denotes a natural number which we take as 3 and 
c is a constant determined with the help of the disc mass. Thus, equation ([2]) becomes: 


1/ = -Gm 


n w ^ _L 


1 Hi 3 H 2 


-Vi 




with ri{i = 1,2) dehne as 


ri = y/(a; + + y"^, r2 = \/(a: + /i - l)^ + y‘^, 


( 3 ) 


( 4 ) 


where G is the gravitational constant. qi and q 2 are the radiation factors of the primaries. 

■p 

Parameter g* = 1 — ^ (i = 1, 2), where iT is the force cause d by radiation pressure force and Eg 


results due to gravitational force fjSingh and Taura l2014all . The oblateness coefficients for the 
bigger primary is denoted as Ai, with 0 < = J 2 iR^ << 1 and for the smaller primary as Bi, 
with 0 < Hj = J 2 iRf‘ « 1, f = 1,2, where J 2 i are zonal harmonic coefficients and Hi ^2 denote 
the mean radii of mi^ 2 , supposing the primaries have their equatorial plai ies coinciding with the 


plane of motion . Vpk is the potential due to the disc which is given by fjjiang and Yeh 


Kushvah et ah 


2 ni 2 h 


2006 


— 4 


E [Q p' {r')r' 


dr', 


( 5 ) 


where E{C,) denotes the elliptic integral of the hrst kind, r' is the disk’s reference radius and 
Q = 2\/rV/(r + r'). We write the Lagrangian to onr problem as 


TTITI / o o\ . . .V TTX / . 0 . o\ - ^ 

L = [x +y )+ mn{xy - xy) +— {x +y ) - Vj^k, 


( 6 ) 
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where the mean motion n is 


— 1 + 2 + Bi — -(^2 + B 2 ) ) - 2/fe(r)|^^^ 


(7) 


fb{r) denotes the gravitational force of the disc fj Jiang and Yeh 120061) 

dVj^k 


fb{r) = 


dr 


0 f pi’f'y 

1 - 

+ 

r=rm Jr' 

/jrt /jrtf rjp 1 ryt' 


dr', 


( 8 ) 


where E{() is elliptic integral of the second kind. We calculate /^(r’) at r = = 0.99. By 

expanding the elliptic integral in equations ([S]) and ([H]) within the limit k < r' < j and c 


an appropriate terms relative to r, we write a simplihed form of /^(r) as (jjiang and Yeh 


loosing 


20061) : 


/b(r) = - 2chTT-^-—^\ - ^chirirLog 

jk r^ 8 H 


(9) 


We assume that the gravitational force is radial 
and y components of the force fbij) respective!} 
mass are: 

X — 2ny = n. 


y symmetric, so we have y/fe(r) and y/6(r) as x 
The equations of motion of the inhnitesimal 


y + 2 nx = kly, 


(10a) 

(lOb) 


where 


o n^{,x^ + y^) , ,, ^ f I ^ 3212^^ , f 1 ^ Bi 3B2\ 

2 s;!)My" ' 

It should be noted here that the suffixes x and y indicate the partial derivatives of hi with respect 
to X and y respectively. This system admits the well-known Jacobi integral: 


c = 2 a-(x^ + f), 


( 12 ) 


where C is the Jacobi constant. 


3 Calculation of libration points 


libration points are very important in astronomy because they indicate places where particle 
can be trapped. In this section, we obtain the libration points which are the singularities of the 
manifold 


C = 2Vt — [x^ + y‘^). 


( 13 ) 
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For detail derivation of equations ©[8] and HI), one is advised to check ref. (|,Tiang and Yeh 1 120061) 
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To achieve this task, we equate all velocities and accelerations of the dynamical systems to zero 
(i.e., X = y = VLx = = ^) ■ Thus, 




dVt 2 ,{3 ~ k\ TlX Sc/lTTX, 

= n X — 2ch 


dx ~ ~ ~ \ jk J 

qi{y - l)(x + y) (-15.42 + 12Air^ + 8rf) 
8r[ 

yq2{x + /i - 1) (-15i?2 + 12i?ir| + 8r^) 

8 rl 


(14a) 




on 


= n^y — 2ch 


j — k\ Try Schiry 


dy y jk J 8r^ 

gi(/i - l)y (-15212 + 122 lir 2 + 8r4) 
^ 8r7 

(-15i?2 + 12i?irf + 8r^) 

8 rl 


Log 


Re-arranging equations fll4ap and I14bp , we hnd 




jk 


8 r'^ 


+ 


yq 2 (-15^2 + l2Birl + 8r|)\ f q,{y - l)y {-IbA^ + + 8rf) 


8 rl 


8r[ 


/i(/i - 1)^2 (-15-02 + 12.Bir| -2 8r|) 
8rl 


= 0 , 


(14b) 


gi(/x - 1) (-15212 + 122liri^ + 8rf) 
8r[ 


(15a) 


y [n^ — 2ch 


j — k\ 71 Sc/iTT 


jk 


8 r^ 


yq2 (-15.02 + 12-0ir| -2 8r2) 

8rl 


Log 


= 0 . 


+ 


gi(/i - 1) (-15212 + 122lir^ + 8r\) 
8r\ 


(15b) 


Now, we solve equations fll5al) and ll5bD for ?/ 7 ^ 0 to obtain the triangular libration points. From 
equation fll5bD . if ?/ 7 ^ 0, then it requires that 


_ 2 ch ilxJA 2 - ^Log 


jk J 8r^ 
yq 2 (-15-02 + 12-0ir^ -2 8r|) 
8 r^ 


+ 


gi(/i - 1) (-I 52 I 2 + 122lir^ + 8rf) 
8r\ 


= 0 . 


(16) 


By substituting equation fITB]) in fll5al] . we hnd 

qi (- 15.42 -2 l2Airl + 8 ri) q 2 (-15.02 + 12 - 0 ir 2 -2 8 r|) 


8r[ 


8 rl 


= 0 . 


(17) 
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Again, we re-write equation ffTHjl in the form of 


— 2ch 


j - k 
jk 


71 

r3 


SchiT 


Log 


gi (-15^2 + 12 Air 2 + 8 rf) 


8 r[ 


-|-/i 


gi (-15^2 + 12 Air 2 + 8 rf) (-15^2 + 125ir2 + 8 r|) 


8 r[ 


8 rl 


= 0 . 


( 18 ) 


Substituting equation ffT7)l into (na, we find 


9 , I i — k\ 71 SchTT^ 

- 2ch i A__ ] _ - —^Log 


jk J 8 r^ 

From equations ffTTD and we obtain 

_ 2ch (^ - ^Log 


jk 


8 r^ 


qi (-15^2 + 12 Airi^ + 8 rf) 

8 rl 


q 2 (-15i?2 + l2Bir2 + 9>r\) 

8 rl 


= 0 . 


(19) 


= 0 . 


( 20 ) 


By neglecting the effects of radiation and oblateness of the primaries, the above equations flT^ 
and fl20|l reduce to classical case with the solutions ri = r 2 = 1. Hence, it is reasonable if we 
assume the solutions of equations flT^ and fl 20 |) as ri = 1 -|- ei and r 2 = 1 -|- e 2 respectively, where 
ei, 62 are very small quantities. Using equations ([7]), ([9]) and flT^ . we hnd 


^ + 77 [ + Bi - -{A 2 + B 2 ) ) - fb{'r)\^ - 


qi (-I 5 A 2 + 12Air2 + 8rt) 


8 rl 


= 0 . 


( 21 ) 


Thus, since e* are very small quantities, hence, (t)/rj, (t)/rf and (t)/rf can all be approximated 
as t (~ t), where t could be any parameter in the numerator of the last expression of equation 
m- Now, substituting gi = 1 — pi, in equation fl2T]l and then restricting ourselves to the linear 
terms in e*, Aj, pi, and coupling terms in qiAi, we have 

Q 1 'i 1 

1 + 2 (^ 1(1 “ '?i) + Bi) —^ (^ 2(1 - qi) + B 2 ) — = ^-Pii (22) 


n = (^l+pi + ^(Ai(l-gi) + Hi)-y (A2(l-gi) + H2)-/b(r)|^^^^^ 

= 1 + “ 2 ^ ^ 8 + B 2 ) + - fb{r)\r^^^ 

= 1 + 61 . (23) 

Using similar approach, 

^2 = 1 + - 2 + ^1) + g (-B2P2 + ^2) + - / 6 (?")lr=r^^ 


1 + 62 . 


(24) 




















Table 1: Effect of the perturbations and PDF on the positions of the triangular libration points 
L 4 5 with /i = 0.4, Tm = 0.99, j = 1.5, k = 1, h = 0.0001 and c = 1910.83 









PDP = 0 

PDP 0 

Case 

Pi 

P2 



A 2 

B 2 

L4,5 

L4,5 

1 

0 

0 

0 

0 

0 

0 

0.100000, ±0.866025 

0.100000, ±0.693135 

2 

0 

0 

0.03 

0 

0.02 

0 

0.102500, ±0.864582 

0.102500, ±0.691691 


0 

0 

0.04 

0 

0.03 

0 

0.101250, ±0.865304 

0.101250, ±0.692413 


0 

0 

0.05 

0 

0.04 

0 

0.100000, ±0.866025 

0.100000, ±0.693135 

3 

0 

0 

0 

0.03 

0 

0.02 

0.108750, ±0.871077 

0.108750, ±0.698187 


0 

0 

0 

0.04 

0 

0.03 

0.110000, ±0.871798 

0.110000, ±0.698908 


0 

0 

0 

0.05 

0 

0.04 

0.111250, ±0.872521 

0.111250, ±0.699630 

4 

0.10 

0 

0 

0 

0 

0 

0.066667, ±0.846780 

0.666667, ±0.673890 


0.15 

0 

0 

0 

0 

0 

0.050000, ±0.837157 

0.500000, ±0.664267 


0.20 

0 

0 

0 

0 

0 

0.033333, ±0.827535 

0.333333, ±0.654645 

5 

0 

0.15 

0 

0 

0 

0 

0.150000, ±0.837157 

0.150000, ±0.664267 


0 

0.25 

0 

0 

0 

0 

0.183333, ±0.817912 

0.183333, ±0.645022 


0 

0.35 

0 

0 

0 

0 

0.216667, ±0.798667 

0.216667, ±0.625777 

6 

0.10 

0.15 

0.03 

0.02 

0.03 

0.02 

0.115417, ±0.821954 

0.115417, ±0.649064 


0.15 

0.25 

0.04 

0.03 

0.04 

0.03 

0.131896, ±0.795071 

0.131896, ±0.622181 


0.20 

0.35 

0.05 

0.04 

0.05 

0.04 

0.148250, ±0.768405 

0.148250, ±0.595514 


Substituting equations fl2l|) and fl2^ into equation (|4]) with restriction to the linear terms in e*, 
we hnd 

1 


Xo = 


yo = 


2 ^ 


1 


- 2/i + 1) = ei - f 

(r^ - (x + + 2 ei - 262)^^ 


1/2 


f( 


4 4 

1 + gf. + 3^2 


1/2 


f( 


2 2 

1 + 361 + 362 


(25) 


Equation (125|) has been used to compute the results in Table [H As it can be seen from the table, 
the coordinates of the triangular libration points L 4 and L 5 are affected by the oblateness and 
the radiation of the primaries when PDP is 0 and 7 ^ 0 . If from equation fl25|) . we neglect coupling 
terms (CT) and PDP (i.e. setting AiPi = BiPi = c = 0, i = 1,2) our results are in excellent 
agreement with the ones obtained in ref. llSiugh and Taura l2014ari for a special case = 0. 
Also, our results coincide with that of ref. fjAbouelmag~ . 2012 1 for Mi, = 0 by neglecting the 
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effects of radiation of the primaries, PDF and oblateness of the less massive primary (i.e. setting 
Pi = 'P‘) = c = Bi = B-) = 0 ). Furthermore, by considering only linear terms in small quantities 
in ref. flKnshvah 1200811 and if we ignore the effects of oblateness of the more massive primary, 
PDP, the radiation and the oblateness coefficient B 2 of the less massive primary in the current 
study (i.e. setting p 2 = A-\ = A 2 . = B 2 = 0 ), the libration points in the current results agree with 
those obtained in ref. (iKnshvah 1200811 for a special case Mb = 0 


For case 1 (as presented in Table [T]), in the absence of the obl ateness, the radia 


tion o 


primaries and PDP, our results agree with that of Singh and Taura (jSingh and Tanra 


the 


201811 . In 


general our results are in excellent agreement with the ones obtained (presented as case 1 — 7) 
by these authors if we set A 2 = B 2 = c = Q except for their cases 2 and 3 which were due to 
erroneous negative in the coefficient of {Ai + A 2 ) in their equation (14). For case 2, we consider 
only the oblateness of the more massive primary (i.e. we set pi = p 2 = Bi = B 2 = 0) and 
in the absence of PDP (i.e. c = 0). The triangular libration points L 4 , L 5 come nearer to 
the primaries and the line joining the primaries. The same behavior is also experience in the 
presence of power prohle. For case 3, we considered only the oblateness of less massive primary 
(i.e. qi = q 2 = l,^i = A 2 = 0). Consequently, the triangular libration points move away from 
the line joining the primaries. 

Generally speaking, the presence of power-law prohle does not affect the x—component of 
the libration points. On considering only the radiation of the less massive primary (i.e. setting 
Pi = Ai = A 2 = Bi = B 2 = 0), the libration points tend towards the line joining the primaries. 
This can be seen in case 4. Furthermore, on considering as case 5, only the radiation of the 
more massive primary (i.e. setting p 2 = Ai = A 2 = Bi = B 2 = 0), the libration points move 
away from the line joining the primaries. As a hnal case, we explored the overall effect due to 
all the perturbations. We hnd that ^ 4 , 5 , move towards the less massive primary. It should be 
noted that cases where power-law prohle is considered exhibit similar behavior with cases where 
power-law prohle is 0 . 

It is worth mentioning that the comprehensive ehects of the perturbations have a stabilizing 
proclivity. However, the oblateness up to J 2 of the primaries and the radiation of the primaries 
have tendency for instability, while coefficients up to J 4 of the primaries have stability predispo¬ 
sition. 


4 Linear Stability 

To study the stability of a libration points (xo,|/o), we apply inhnitesimal displacement C and rj 
to the coordinates by using p = y — Po and ( = x — Xq and then substitute them into equations 
of motion in fllOal) and (llObD to obtain 


C- 2 np = CDL + 


xy 


p + 2nc = cnl^ + pa 


yy' 


(26) 
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where the superhx ‘O' indicates that the partial derivatives have been computed at the triangular 
libration points by considering {xo,yo). Now, let us assume a solution of the form ( = Ciexp"^' 
and r] = C 2 exp^*, where Ci and C 2 are constants and A is a parameter. Substituting the assumed 
solutions into equation fl26|) . we obtain the following non-trivial solutions for Ci and C 2 


y - fiL -2nA - ai, 

2 n\ - 


0 . 


(27) 


Solving the determinant by expansion, the characteristic equation corresponding to the varia¬ 
tional equations fl26|) can be found as: 


A- + (4n^ - jiL - + (nlX, - = 0 . 


(28) 


Now, we hnd the second order derivatives of the potential function w.r.t x and y and we write 
the expressions as follows: 
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To get and we substitute the libration points given by equations fl25|l into equations 

fl29a[ I29bl and I29cp . Furthermore, we restrict our calculations to only the linear terms in e*, Ai, 
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Pi and coupling terms in qiAi and qiBi. Thus, we obtain 
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More explicitly, we can re-write the above expressions for VPyy and 17°^ as 
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Solving the characteristic equation fl25]) for A, we obtain 
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where the discriminant 6 is given by 
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It can be deduced from fl36|) that (5|^^g > 0 and (3|^^i < 0 (i.e., the values of 6 at the point 
p = 0 and /i = I are of opposite signs). Also from equations fl371) and fl38|l . <5^ < 0 within the 
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interval of (O, |). Since S is strictly decreasing fnnction of jj, in the open interval of (O, |), thns, 
there is only one valne of /i in the interval for which 6 vanishes and it is denoted as jUc (the 
critical mass parameter). Moreover, there exist three possible cases. When 0 < < jUc (3 > 0), 

the valnes of Ai, 2 , 3,4 which are the solntions of the characteristic eqnation are pnre imaginary 
nnmbers indicating that, triangnlar equilibrinm points are linearly stable. Secondly, for /i = /ic 
(5 = 0), then Ai ,3 = A 2,4 (i.e., donble roots). The solntions have a secnlar terms which gives the 
instability of the point. Finally, for /Xc < /x < ^ (5 < 0), Ai, 2 , 3,4 have imaginary roots and two of 
the roots have positive real parts. Conseqnently, the triangnlar points are nnstable. In the next 
section, we proceed to calculating pc. 


5 Determination of 


A critical mass dehne a range for which the triangular libration points are stable or unstable. 
It could be a lower or upper limit for this range of values. To calculate the critical mass, we 
re-write the discriminant 5 in the form 
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It is very difficult to obtain an analyt ical solution of equation ([39]). To overcome this diffi¬ 
culty, Abouelmagd (lAbonelmagd 120131) developed an a l gorith m to perform the calculation in a 
simplihed form. Therefore, following ref. ( Abouelmagd . 20131) . we calculate /ic as follows: 


- 


27(1 + - V729(1 + Cm,)^ - 108(1 + C^J(1 + C^ 3 ) 


54(1+ 


-(1-hC -C ) - il - C ) 

' '^rn2 '-'mi) '-'mi) 
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^ 621 ~ 6 ^ 


2 (^ ■ ~j + -C™.) + 


127C„, - 2C„3 - 26C„,| 


(43) 


The expression | (^1 ~ ~ 0.0385... denotes Routh criterion. It determines the stability 

condition for in classical R3BP. The remaining expressions in the Hc equation are due to the 
effects from the oblateness of the more massive and less massive primaries, radiation pressure 
of the less and more massive primaries and PDF. The critical mass (in 17 decimal places) can 
be computed via the substitution of htting parameters = 0.99, j = 1.5, k = 1, h = 0.00001, 
c = 1910.83. Consequently, we hnd 


jjc = 0.03086268515886475 - 0.09804047314470574Ai + 0.8293666279034516Aipi 

-0.08625060012103034Ri - 4.17143645655834750RiPi + 0.36284935058668233R2P2 
-0.14427151017244513pi + 0.39017449819889993^2 - 0.98689512656254470^2^1 
-0.009527626111967606^2 + 0.19230243478937234R2 + 5.23865406035375200R2Pi 
-0.3230170771089508 AiP 2 + 0.41654634990801520^2^2 + 0.02473745338809197Rip2. 

(44) 


We now proceed to numerically compute the criti cal mass for th e Eart h-Moon and Jupiter-Moons 
systems by using the astrophysical parameters (lAbonelmagd I . l2012l) presented in Table [2] with 
Ai = Bi and A 2 = - 82 - Our results are displayed in Tabled The PDP has effect of about ~ 0.01 
reduction on the application of /ic to Earth-Moon and Jupiter-Moons systems. The result reflects 
how ring of dust which could be as a result of asteroid collisions can affect the stability range 
of an inhnitesimal body in R3 Bf|^. Furthermore, in hgure [21 we have plotted the variation of 
the critical mass ratio as a function of A 2 and Ai and then varies other parameters in each case. 
The zonal harmonic coefficients decrease and increase the critical mass in respect to a classical 
case (i.e. when Ri = A 2 = i ?2 = Pi = P 2 = 0). As it is anticipated, in both hgures (a and b), it 
is observe that the critical mass is less than 0.084. 
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For instance, stability range of a spacecraft in Jupiter-Callisto system. 
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Table 2 : Astrophysical parameters for planets systems of our solar system 


S.N. 

System 

Ai 

A 2 

Earth 

1 

Moon 

0.000000298147 

0.000000000000 

Jupiter 

2 

Callisto 

0.000021186102 

-0.000000001213 

3 

Metis 

0.004586424842 

-0.000056862695 

4 

Adrastea 

0.004515511583 

-0.000055117917 

5 

Pasiphae 

0.000000136024 

0.000000000000 

6 

Sinope 

0.000000133738 

0.000000000000 


6 Discussions and conclusions 


The results obtain in this study can be reduce to the ones obtained previously if an adjustment 
is made to the model parameters by setting some of them to zero. For instance, the equations 
of motion are in agreement with 


Classical case of Szebehely fjSzebehelv 1.119671) . if the effects of perturbations and power-law 
density prohle of disc are ignored (i.e. setting = p 2 = = ^2 = -Bi = -^2 = c = 0 ) 


The result of Singh and Taura (jSingh and Taura I . l2013l) for Mf, = 0 if the effects of power- 
law density prohle of disc (PDF) are ignored (i.e. setting c = 0). 


The result of Singh and Taura (ISingh and Taura I . I2ni4bfl for Mb = 0 in the absence of the 
oblateness up to J 4 of the primaries and PDP (i.e. setting A 2 = B 2 = c = 0). 


The result of Kishor and Kushvah flKnshvah et ah 1 . 120121) . if we neglect the effects of the 
radiation factor of the less massive primary, oblateness coefhcients of the more massive 
primaries and up to J 4 of less massive primary (i.e. setting Ai = A 2 = B 2 = q 2 = 0). 


Abouelmagd flAbonelmagd 1 . 1201211 . if we consider the oblateness coefhcients up to J 4 and 
absence of PDP of the more massive primary (i.e. setting pi = p 2 = Bi = B 2 = 0), 


Singh and Ishwar fISingh and Ishwar llOOOh . in the presence of the oblateness up to J 2 of 
the primaries and in absence of PDP (i.e. setting A 2 = B 2 = c = 0). 
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Table 3: Critical mass for Earth-Moon and Jupiter-Moons system. 


Power-law prohle = 0 

Power-law prohle 7 ^ 0 

Pc (Earth-Moon) 

Pc(Jupiter-Calhsto) 

Pc(Earth-Moon) 

Pc(Sun-Calhsto) 

0.03852079631 

0.03851353085 

0.03086264505 

0.03080640048 

Pc(Jupiter-Metis) 

Pc( Jupiter-Adrastea) 

Pc(Jupiter-Metis) 

Pc (Jupiter-Adrastea) 

0.03687828788 

0.03690440886 

0.02998434157 

0.02994725276 

Pc(Jupiter-Sinope) 

Pc (J upiter-Pasiphae) 

Pc(Jupiter-Sinope) 

Pc(Jupiter-Pasiphae) 

0.03852085349 

0.03852085269 

0.03081027497 

0.03086267535 


Kushvah flKushvah I . l2008l) when only linear terms in small quantities are considered with 
the effect of the radiation of the more massive primary, oblateness up to J 2 of the less 
massive primary are considered and ignoring the effect of PDF (i.e. setting q 2 = l,Ai = 
A 2 = B 2 = c = 0). 


Singh and Taura (jSingh and Taura I . l2014bl) . if we disregard the effect of the radiation of 
the primaries and PDP (i.e., setting Pi = P 2 = c = 0). 


Jiang and Yeh fj Jiang and Yeh I . l2006l) in the absence PDP and all the perturbations (i.e. 
setting = p 2 = Ai = A 2 = Bi = B 2 = 0). 


In the limiting case where control factor of the density prohle (c) equals/tend to zero, then, pc 
obtained in equation fl44|) becomes 


Pc = 0.03852089650455137 - 0.285001787790555607li + 0.69578421210336550712 

-0.062779565568333305i + 0.1402286565478101452 + 0.3845227397304622552^2 
-0.00891747059894587pi + 0.629660595060427707lipi - 0.49510688989916757l2Pi 
-3.756468386787854005ipi + 4.7178791599821820052Pi - 0.00891747059894587^2 
-0.312023942343409607liP2 + 0.41232360442662697712^2 - 0.004358682110150005iP2 

(45) 


It is worth mentioning that if we neglect coupling terms AiPi and BiPi {i = 1, 2) in the above 


calcula tion, the result reduces to the one obtained previously in Eq. 129 1 of ref (ISingh and Taura 


2014a 1 for = 0 and also in excellent agreement with Eq. (21) of ref Jsingh and Taura . 2014b 1 
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Figure 2: (Color online) (a) Critical mass ratio as a function of A 2 with (I) Ai = 0.03, Bi = 0.02, B 2 = 0.01, 

<71 = 0.99 and <72 = 0.98. (II) Ai = 0.03, Bi = 0.02, B 2 = 0, qi = 0.99 and (72 = 0.98. (Ill) Ai = Bi = B 2 = 0, 

<71 = 0.99 and (72 = 0.98. (IV) Ai = Bi = B 2 = pi = P 2 = 0. (b) Critical mass ratio as a function of Ai with 

(I)i3i = 0.02 A 2 = 0.02, B 2 = 0.01, qi = 0.99 and (72 = 0.98. (II) Bi = 0.02 A 2 = B 2 = 0, <71 = 0.99 and 

<72 = 0.98. (Ill) Bi = A 2 = B 2 = 0, <71 = 0.99 and (72 = 0.98. (IV) Bi = A 2 = B 2 = pi = P 2 = 0. 


for Mh = 0. Furthermore, if we neglect the effect of radiation of the primaries, CT, and oblateness 
of the less massive primary (i.e. setting pi = V 2 = i?i = i ?2 = 0 ), then equation (HSj) reduces 
to the one obtained in ref. flAbouelmagd I . l 2012 l ). Again, if only the oblateness up to J 2 of the 
primaries is conside red and ignoring CT (i.e. A2 = B2 = 0), then eqnation 0451) rednced to the 
one obtained in ref. fISingh and Ishwar 1 . 119991) . Moreover, by ignoring effect of all pertnrbations 
in eqnation (H5D (i.e. setting pi = p? = A^ = Ao = = i ?2 = 0 ), then onr resnlt becomes that 

of classical CR3BP presented in ref fjSzebehelv 119671) 

In this paper, we have stndied the effect of oblateness np to J 4 of the primaries within the 
framework of restricted three body problem, by nsing a more realistic model in which a disc 
with PDF is rotating aronnd the common center of the system mass. The existence and stability 
of triangnlar eqnilibrinm points have been explored. It has been demonstrated that triangnlar 
eqnilibrinm points are stable for 0 < p < p-c and nnstable for Pc < P < 1 / 2 , where pc denote 
the critical mass parameter. Conseqnently, the oblateness np to J 2 of the primaries and the 
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radiation can rednce the stability range, while the potential from the circnlar clnster of material 
points and oblateness np to J 4 of the primaries will increase the size of stability both in the 
context where PDF is considered and ignored. The PDP have effect of abont 0.01 rednction on 
the application of fic to the Earth-Moon and Jnpiter-Moons system. In the limiting case c = 0, 
and also by setting appropriate parameter(s) to zero, onr results are in excellent agreement with 
those obtained previously. Our result have a practical application in sp ace dynamics. 

Furthermore, the results obtained in ref. (lAbonelmaEfd et ah 120151) have been extended and 
generalized in the current study, by considering the effect of PDP and radiation and radiating 
factor on the triangular libration points and the critical mass within the framework of R3BP. 
Our results reflect how ring of dust which could be as a result of asteroid collisions can affect 
the stability range of an inhnitesimal body in R3BP. 

All our numerical computations can be done using a commercial symbolic package such as 
Mathematica, Maple or MATLAB etc. 

Finally, we suggest a possible extension of the current work to a case of non-coplanar equi¬ 
librium. 


Acknowledgments 

We thank the kind referees for the positive enlightening comments and suggestions, which have 
greatly helped us in making improvements to this paper. In addition, BJF acknowledges eJDS 
(ICTP). OAF acknowledges University of Ilorin for granting him sabbatical leave. The Astro¬ 
physics group acknowledges supports from the Vice Chancellor of Federal University Laha, Prof. 
Ekanem Ikpi Braide, Dean of Science Faculty, Prof. Emmanuel Kwan-Ndong and the Head 
of Physics Department, Prof. Sanusi Mohammed Liman. This work is supported partially by 
20150964-SIP-IPN, Mexico. 


References 

Abouelmagd, E. I. Existence and stability of triangular points in the restricted three-body prob¬ 
lem with numerical applications. Astrophys. Space Sci. 342 ( 2012 ) 45. 

Abouelmagd, E. I. Stability of the triangular points under combined effects of radiation and 
oblateness in the restricted three-body problem. Earth Moon Planets 110 (2013) 143. 

Abouelmagd, E. L, Awad, M. E., Elzayat, E, M. A., Abbas, I. A. Reduction the secular solution 
to periodic solution in the generalized restricted three-body problem. Astrophys. Space Sci. 
350 (2014) 4 


22 





Abouelmagd, E. L, Alhothuali, M. S., Guirao, J. L. G., Malaikah, H. M. The effect of zonal 
harmonic coefficients in the framework of the restricted three-body problem. Adv. Space Res. 
55 (2015) 1660. 

Abouelmagd, E. L, Mostafa, A. Out of plane equilibrium points locations and the forbidden 
movement regions in the restricted three-body problem with variable mass. Astrophys. Space 
Sci. 357 (2015) 58. 

Abouelmagd, E. I., Sharaf, M. A. The motion around the libration points in the restricted three- 
body problem with the effect of radiation and oblateness. Astrophys. Space Sci. 344 (2013) 
321. 

Falaye, B. J. Effect of oblateness, radiation and a circular cluster of material points on the 
stability of equilibrium points in the restricted four-body problem. Few-Body Syst. 56 (2015) 
29. 

Jiang, I. G., Yeh, L. G. Galaxies with supermassive binary black holes: (I) A possible model for 
the centers of core galaxies. Astrophys. Space Sci. 349 (2014) 881. 

Jiang, I. G., Yeh, L. C. On the Ghermnykh-like problems: I. The mass parameter p = 0.5. 
Astrophys. Space Sci. 305 (2006) 341. 

Kholostova, O. Stability of triangular libration points in a planar restricted elliptic three body 
problem in cases of double resonances. Int. J. Non-Linear Mech. 73 (2015) 64. 

Kushvah, B. S. Linear stability of equilibrium points in the generalized photogravitational 
Chermnykh’s problem. Astrophys. Space Sci. 318 (2008) 41. 

Kushvah, B. S., Kishor, R., Dolas, U. Existence of equilibrium points and their linear stability in 
the generalized photogravitational Chermnykh-like problem with power-law profile. Astrophys. 
Space Sci. 337 (2012) 115. 

Letelier, P. S. da Silva, T. A. Solutions to the restricted three-body problem with variable mass. 
Astrophys. Space Sci. 332 (2011) 325. 

Narayan, A., Shrivastava, A. Pulsating different curves of zero velocity around triangu¬ 
lar equilibrium points in elliptical restricted three-body problem. J. Math. (2013) doi: 
10.1155/2013/936859. 

Singh J., Cyril-Okeme, V. U. Perturbed robe’s circular restricted three-body problem under an 
oblate primary. New Astronomy 34 (2015) 114. 


23 


Singh, J., Ishwar, B. Stability of triangular points in the generalised photogravitational restricted 
three body problem. Bull. Astron. Soc. India 27 (1999) 415. 

Singh, J., Taura, J. J. Combined effect of oblateness, radiation and a circular cluster of material 
points on the stability of triangular libration points in the R3BP. Astrophys. Space Sci. 351 
(2014) 499. 

Singh, J., Taura, J. J. Effects of zonal harmonics and a circular cluster of material points on the 
stability of triangular equilibrium points in the R3BP. Astrophys. Space Sci. 350 (2014) 127. 

Singh, J., Taura, J. J. Motion in the generalized restricted three-body problem. Astrophys. Space 
Sci. 343 (2013) 95. 

Szebehely, V. : Theory of orbits: The restricted problem of three bodies. Academic Press, San 
Diego (1967). 

Trivailo, O. Spacecraft stability, dynamics and control near the triangular Lagrange points in¬ 
fluenced by multiple Trojan asteroids. Acta Astron. 60 (2007) 658. 


24 


